Report for Model 4

— for both jointly and separately

2024-08-10

1. Joint Modeling

1.1 Trace plots for convergence check

The current MCMC setting is:

  • 100,000 iteration;
  • 90,000 burn-in;
  • 10 thinning.

1.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.00       1.00
## LDevsum       1.01       1.03
## dh0           1.25       1.68
## dh1           1.25       1.68
## dh2           1.00       1.01
## dl0           1.06       1.18
## dl1           1.18       1.47
## 
## Multivariate psrf
## 
## 1.18

1.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

1.4 WAIC results

LevelH LevelL
DIC 1141.77712 22552.0905
DIC3 1111.55453 22638.2755
PWAIC 22.66121 203.6814
WAIC 1122.82311 22651.9122

2. Separate Modeling

2.1 Trace plots for convergence check

The current MCMC setting is:

  • 100,000 iteration;
  • 90,000 burn-in;
  • 10 thinning.

2.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.01       1.04
## LDevsum       1.00       1.01
## dh0           5.56      11.07
## dh1          10.17      18.41
## dh2           5.84      10.29
## dl0           1.11       1.32
## dl1           1.02       1.05
## 
## Multivariate psrf
## 
## 9.02

2.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

2.4 WAIC results

LevelH LevelL
DIC 1388.1021 33774.6195
DIC3 1244.4182 33891.6284
PWAIC 89.7843 169.8478
WAIC 1312.5302 33896.8859